More Reliable Measurements of the Slip Length with the Atomic Force Microscope 
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Further improvements are made to the non-linear data analysis algorithm for the atomic force 
microscope [P. Attard, arXiv: 1212.3019 v2 (2012)]. The algorithm is required when there is curva- 
ture in the compliance region due to photo-diode non-linearity. Results are obtained for the hy- 
drodynamic drainage force, for three surfaces: hydrophilic silica (symmetric, Si-Si), hydrophobic 
dichlorodimethylsilane (symmetric, DCDMS-DCDMS), and hydrophobic octadecyltrichlorosilane 
(asymmetric, Si-OTS). The drainage force was measured in the viscous liquid di-n-octylphthalate. 
The slip- lengths are found to be 3nm for Si, 2nm for DCDMS, and 2nm for OTS, with an uncer- 
tainty on the order of a nanometer. These slip lengths are a factor of 4-15 times smaller than those 
obtained from previous analysis of the same raw data [L. Zhu et al., Langmuir, 27, 6712 (2011). 
Ibid, 28, 7768 (2012)]. 



I. INTRODUCTION 

Measurements of the slip length are significant for 
three reasons. First, the hydrodynamic equations for any 
flow can't be solved without specifying the boundary con- 
ditions, and so whether a fluid sticks or slips during shear 
flow at a solid surface is fundamental to the application of 
hydrodynamics to the real world. In turn, the extent that 
a fluid slips, if any, effects quantitatively a range of phys- 
ical phenomena (e.g. flow rates in pores, drainage forces 
between particles, lubrication of surfaces, pressure heads 
in microfluidic devices), and so measuring, understand- 
ing, and controlling slip could lead to new technologies, 
devices, and industrial processes. 

Second, whereas equilibrium statistical mechanics is 
well-established for elucidating static molecular structure 
(in bulk and at surfaces), the same cannot be said for 
the non-equilibrium case, either in general or for fluid 
flow. Hence nanoscopic measurements of the slip length 
yield fundamental information about the structure of the 
fluid and its interaction with the solid surface during 
shear flow. Such data is valuable in its own right, pro- 
viding molecular-level insight into inhomogeneous flow 
and perhaps enabling one to identify the specifically non- 
equilibrium aspects of the way in which structure, flow, 
and interaction are entwined at the molecular level. In 
turn the measurements provide specific motivation to de- 
velop non-equilibrium computer simulations, and bench- 
marks against which those algorithms could be tested. 

Third, slip can be regarded as a correction to stick 
boundary conditions, and as such it is a second order ef- 
fect that is a real challenge to measure experimentally 
with any reliability or accuracy. Obviously any small 
error in the measurement overall translates into a large 
error in the slip length. The literature abounds with 
examples where the reported slip lengths vary by an or- 
der of magnitude for ostensibly the same system. ^'^ (The 
present paper will present slip lengths that are a factor of 
15 smaller than those previously reported for the same 
data.) Hence credible measurements of the slip length 
can only be obtained by improving the reliability and 
accuracy of the measurement technique itself, and this 
would have broad benefits beyond the immediate appli- 



cation to slip. 

On this last point, the present paper uses the atomic 
force microscope to measure the slip length via the anal- 
ysis of the hydrodynamic drainage force. The present 
results represent the third generation of improvements 
to atomic force microscopy that have been motivated by 
the desire for reliable measurements of the slip length. 
Most of the new algorithms and procedures are general 
in nature and can be directly applied to other forces mea- 
sured with the atomic force microscope. In the present 
case, where the sought for slip length is a second order ef- 
fect, the new procedures have proved essential for reliable 
results. In other cases, where one might be measuring a 
first order effect, the improvements are still worthwhile 
because they reduce the quantitative error of the mea- 
surements by about an order of magnitude while still re- 
maining simple enough for the data analysis to be carried 
out by a spreadsheet. 

In Appendix]^ the various improvements in force mea- 
surement and data analysis with the atomic force micro- 
scope are summarized. The focus is on those develop- 
ments stimulated by the challenge of accurate and reli- 
able measurements of the slip length, although of course 
many of these improvements are more generally appli- 
cable. As discussed in the appendix, the most recent 
series of improvements concern the non-linear analysis 
of the measured force data^i The need for such an algo- 
rithm arises when the compliance region has curvature, 
as is shown in Fig. [T] The non-linearity is due to the re- 
sponse of the photo-diode, most likely arising from non- 
uniformity in intensity and width of the light beam mov- 
ing across the split photo-diode. In Ref. |j, an algorithm 
was given based on a polynomial of best fit to the contact 
region, from which the angular and vertical deflection of 
the cantilever in the non-contact region can be obtained. 
The algorithm was applied to measured raw data for the 
drainage force in the Si-DOPC-Si system. 

This paper presents further improvements to the non- 
linear data analysis algorithm that are applicable to gen- 
eral atomic force microscopy. The modifled analysis is 
applied to two new systems, DCDMS-DOPC-DCDMS 
and Si-DOPC-OTS, as well as re-analysing the Si-DOPC- 
Si system. 
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FIG. 1: The raw photo-diode voltage versus the piezo-drive 
displacement. The solid curve is measured extension data, the 
dashed line gives the tangent at first contact, and the dotted 
line gives the average slope in contact. The inset shows the 
analysed force versus separation, with the solid curve resulting 
from the non-linear analysis and the dashed and dotted curves 
resulting from the conventional linear analysis using the first 
contact slope and the average slope, respectively. The source 
of the measured data is Ref. 9, which may be consulted for 
the experimental details. 

One improvement concerns replacing an extrapolation 
of the polynomial fit by an inherently more reliable inter- 
polation. This can be explained as follows. The voltage 
measured on extension can be split into two ranges, in 
contact V > T4,oxt, and prior to contact V < Vc,ext (be- 
cause the drainage force is repulsive), where T4,ext is the 
voltage at initial contact. This means that in order to 
obtain the pre-contact drainage force on extension from 
the measured voltage, one has to use an extrapolation 
of the non- linear fit made in contact. This extrapolation 
can lead to small but unacceptable errors. In contrast, 
on retraction, in contact the voltage is > Vc,rct, and 
after contact V ^ V^,rot (because the drainage force is 
attractive), the voltage at final contact, Vc^rct, is strictly 
less than the voltage measured on the whole range of 
extension. This means that one can use the non-linear 
fit to the retraction data in contact to obtain both the 
extension and the retraction out of contact force by in- 
terpolation rather than extrapolation. This gives a small 
but significant improvement in the drainage results. 

A second improvement is that a least squares fit algo- 
rithm has been developed to obtain the cantilever spring 
constant and effective drag length from the measured 
data at large separations. By automating this proce- 
dure, a source of possible human bias is eliminated and 
the time required to analyze each force curve is signifi- 
cantly reduced. 

A third improvement was also explored, which in some 
systems can be important, but which makes negligible 
difference for the atomic force microscope data analyzed 
here. This is variable cantilever drag, which arises from 
the fact that the hydrodynamic drag on the cantilever de- 
creases as the cantilever bends in response to the drainage 



force. This effect was not taken into account in the theo- 
retical forces that were fitted to the non-linearly analyzed 
dataj^ In earlier linear analysis, variable drag was found 
to be significant for weak cantilevers {kcs ^ .1 N/m), but 
negligible for stiff ones(fceff ^ 1 N/m).^"~ Since the spring 
constant obtained with the non-linear analysis was rela- 
tively large, fccff = 1.68 N/m,"* it was assumed that it 
would be acceptable in the first instance to neglect vari- 
able drag. In the present work this assumption is checked 
by performing variable drag calculations for the experi- 
mental conditions, and it is found that it does indeed 
have negligible effect. However in the process of per- 
forming the check, a more robust numerical algorithm 
was developed that gives more reliable results and that 
allows the retract data to be calculated as well. This 
algorithm has some intrinsic interest and is included in 
Appendix B, which sets out in detail the non- linear data 
analysis algorithm. 

For an independent experimental test of the present 
protocols, the procedure used in Ref. H is used here 
as well. In that case three sets of measurements were 
performed for the drainage force in the liquid DOPC 
(di-n-octylphthalate): Si-Si, DCDMS-DCDMS, and Si- 
OTS, where Si is a silicon wafer with a native silicon 
oxide layer, DCDMS is a dichlorodimethylsilane self- 
assembled monolayer prepared on the same type of sili- 
con wafer from the vapor phase, and OTS is an octade- 
cyltrichlorosilane self-assembled monolayer also prepared 
on a silicon wafer. In water. Si is hydrophilic with a con- 
tact angle close to zero, whereas DCDMS is hydrophobic 
with an advancing contact angle of 109°,^ and OTS is also 
hydrophobic with contact angle of 112°.^ In DOPC the 
contact angle is 21° for Si^^ 48° for DCDMS,^ and 45° 
for OTSj^ Given the similarity of the two hydrophobic 
monolayers, one would expect them to have similar slip 
lengths, and the two slip lengths fitted for the symmet- 
ric systems ought also fit the asymmetric system. This 
represents a test of the reliability of the measurement 
protocol. 



II. RESULTS 

Raw data for three systems were analysed: Si-Si, 
DCDMS-DCDMS, and Si-OTS. These are summarised in 
Table m The tilt angle of the cantilever was -11°. The 
cantilever length quoted in the table is Lq — Lq — 2R, 
where Lq is the nominal length quoted by the manufac- 
turer. This recipe accounts approximately for the fact 
that the colloid probe is mounted slightly back from the 
leading edge of the cantilever. 

There is a one-to-one relationship between the intrinsic 
cantilever spring constant kg and the effective force mea- 
suring spring constant kes, which is given in Eq. (jC33[) 
below. Usually both are quoted below even though this 
is redundant. In brief, the effective spring constant gives 
the surface force from the vertical deflection of the con- 
tact point, Fz = koSZc, whereas the intrinsic cantilever 
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TABLE I: Parameters and results for three series of drainage 
force measurements.* 





Si-Si 


DCDMS-DCDMS 


Si-OTS 


No. Meas. 


9 


11 


10 


Cantilever^ 


C 


A 


F 


Lq (/xm) 


110 


90 


230 


R (/im) 


10.11 


9.28 


10.14 


fco (N/m) 


1.38 ±.06 


1.78 ±.10 


1.20 ± .07 


fceff (N/m) 


1.69 ±.08 


2.22 ± .12 


1.34 ± .08 


-i/drag (Atm) 


85 ±8 


71 ±5 


162 ±5 


A (Jm-2) 


1 X 10~^^ 


5 X lO"^" 


1 X 10"^° 


6 (nm) 


3 


2 


3,2 



*A11 measurements were performed in di-n-octylphthalate 
(DOPC), with viscosity -q = 50-54 mPas. The average and 
standard deviation over different drive velocities are given. 
Here Lq is the length of the cantilever, R is the radius of the 
colloid probe, the k are cantilever spring constants (see 
text), Ldrag is the effective drag length, A is the Hamaker 
constant, and h is the slip length. 
^NSC12, tipless, rectangular (Mikromasch). 
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FIG. 2: Vertical cantilever deflection Zc versus separation h at 
a drive rate of ip = 50/ims~^ for Si-Si (see Table The ef- 
fective cantilever spring constant is k^a = 1.69 N/m. The cir- 
cles are the analysed experimental measurement (every tenth 
point plotted), the dotted curve is the calculation with stick 
boundary conditions, b = Onm, and the dashed curve is the 
slip calculation, b = 3 nm. The source of the raw measured 
data is Ref. [3. 



spring constant gives the hypothetical force acting nor- 
mal to the axis of the cantilever that would cause a deflec- 
tion of the cantilever also normal to its axis, F±_ = ko(. 
Most calibration techniques yield fco and this is what is 
usually quoted and used in papers even though it is fccff 
that is required for the quantitative analysis of atomic 
force microscopy. In the present paper the spring con- 
stant was determined by a least squares fit of the drainage 
and drag force to the measured data at large separations. 
The error is about one third of the standard deviation 
shown in Table ID 

In the figures below the vertical deflection of the con- 
tact point Zc is generally given. In previous papers this 
was called the vertical deflection of the tip and was de- 
noted Zt^'^ It is more precise to call it the vertical de- 
flection of the contact point, and to instead write ClCq 
for the vertical component of the deflection of the end 
of the cantilever (see 31) . In contact, Azc = —Azp, 
where Zp is the piezo-drive position. As just mentioned, 
the surface force is related to this by Fz = kcsZc. 

In the deflection-separation curves given below, the 
constant cantilever drag force has been subtracted. The 
theoretical calculations also have constant drag sub- 
tracted. Calculations have been made with constant 
and with variable cantilever drag. For the three systems 
treated below, the effect of variable drag is negligible. 
The variable drag algorithm used in places below im- 
proves upon the original^ in that it is more robust (in 
fact, it is more reliable than the spread sheet algorithm 
for constant drag), it includes tilt and torque, and results 
are obtained for adhesion and retraction. 



A. Si-Si 

Figure [2] shows atomic force microscope measurements 
of the drainage force for Si-Si. The vertical cantilever de- 
flection is shown, from which the force may be obtained 
as Fz — kcsZc- The raw data was analysed using the 
non-linear procedures described in Appendix |B] It can 
be seen that the drainage force is repulsive on exten- 
sion (approach) and attractive on retraction (retreat). 
At large separations there is almost complete overlap be- 
tween the measured data and the two calculated curves. 
At small separations, it can be seen that the stick the- 
ory significantly overestimates the magnitude of the mea- 
sured drainage force on both extension and retraction. 
The slip theory with 6 = 3 nm shows an almost perfect 
fit on extension all the way into contact. The agreement 
is less good on retraction immediately pulling out of con- 
tact, but by about 40 nm the measured and calculated 
deflection converge. 

The same raw data was non- linearly analysed in Ref. 0- 
The difference in the two is that the present analysis uses 
the non-linear flt for the voltage in contact on retraction 
to obtain the cantilever deflection angle out of contact 
on extension. Also, here the cantilever spring constant 
was obtained by a least squares flt of stick theory to 
the large separation data, h G [0.9, 3.7] /xm. This was 
done at each drive velocity and the result averaged over 
all velocities. In Ref. the flt was performed by eye. 
The value of the spring constant obtained previously was 
kcff = 1.68 N/m;^ whereas the present least squares flt 
gave fcoff = 1-69 ± .08 N/m, which suggests that the au- 
thor has a very good eye. The slip length obtained here, 
6 = 3 nm, is the same as that obtained previously.^ In 
the original linear analysis of the atomic force microscope 
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FIG. 3: Deflection Zc at first contact versus drive velocity 
for Si-Si. The symbols are the analysed measured data (raw 
data from Ref. |3|), the dotted curve is the stick calculation 
h = Onm, and the dashed curve is the slip calculation with 
h = 3nm, with other parameters as in Table [H The bars are 
a crude estimate of the error arising from the uncertainty in 
identifying first contact. 



data, the spring constant fitted was fcos = 1.5 N/m, and 
the slip length fitted at low shear rates was &o = 10 nm."^ 

The slip theory fits the measured data on extension 
almost perfectly in Fig. [21 and this level of agreement 
is reasonably typical of all the forces analysed for Si-Si. 
This can be seen in Fig. [3l where the cantilever deflection 
at first contact on extension is shown. In the case of the 
experiments, there can be an uncertainty (on the order 
of several nanometers in the drive distance) in deciding 
exactly where first contact occurs, and it is crudely esti- 
mated that this gives an error on the order of 2-10 nm in 
the deflection. Alternatively, the lack of smoothness in 
the measured data points in Fig. [3] also gives a guide to 
the error in the deflection at flrst contact. Within this 
error, it can be seen that a slip length of & = 3 nm gives a 
deflection at first contact in quantitative agreement with 
the measured one over the whole range of drive velocities. 
The stick theory results shown in the figure correspond 
to h ~ Onm and overestimate the deflection by about 
35% at z — —50 /im/s. From this one can conclude that 
any variation in the slip length by more than a fraction 
of a nanometer will lead to a signiflcantly worse fit than 
has been obtained for 6 = 3 nm. 

The van der Waals force has very little effect on the 
force on extension. For example, at a drive velocity of 
z = — 50/zm/s, the calculated defiection at first con- 
tact for 6 = 3 nm and a Hamaker constant oi A = 
1.3 X 10^^^ J m^^ is 101. 4nm. Changing the Hamaker 
constant to A = 5 x 10^^^ J m^^, the deflection becomes 
100.9 nm. 

Figure S] shows the adhesion deflection, which is de- 
fined as the negative of the minimum deflection that oc- 
curs on retraction. This is easier to obtain and shows 
less variability than the deflection at last contact. In 
general terms obtaining the adhesion reliably in atomic 



FIG. 4: Adhesion defiection Zc versus drive velocity for Si-Si. 
The symbols, curves and parameters are as in the preceding 
figure. The adhesive tension is = kcaZc- The calculations 
use a Hamaker constant of A = 1.3 x 10^^^ Jm~^ and zero 
force position of zo = 0.53 nm. 

force microscopy is extremely challenging. This is be- 
cause the jump out of contact is a catastrophic event and 
it is exceedingly sensitive to external vibrations, surface 
roughness, elastic deformation, and to sliding and rolling 
of the probe, and to peeling of the contact region. The 
present drainage force appears to be an exception, since 
there is quite good agreement between theory and mea- 
surement in this case. Although the instant of the jump 
out is variable, the minimum in the force curve, which is 
the maximum tension, is quite stable and appears to be 
largely determined by the drainage force. 

It can be seen in Fig. S] that stick theory significantly 
overestimates the adhesion, whereas slip theory with b = 
3 nm is relatively accurate. From the high velocity data 
one might conclude that 6 = 3 nm is an upper bound on 
the slip length. 

Although the actual pull-off deflection at last contact 
can be sensitive to the van der Waals force, the adhesion 
deflection as defined here is less so. For example, at a 
drive velocity of i = —50 /im/s, the measured adhesion 
deflection is 98.3 nm. The calculated adhesion deflection 
for a Hamaker constant of A = 1.3 x 10~^"'^Jm~^ and 
a slip length of 6 = 3nm is 86.7 nm, and that for stick, 
6 = Onm is 136.0 nm. Changing the Hamaker constant 
to A — 5 X 10~^^ J m~^ for the slip length of 6 = 3 nm, 
the deflection becomes 88.9 nm. 



B. DCDMS-DCDMS 

Analysed atomic force microscope results for the verti- 
cal defiection of the cantilever as a function of separation 
for DCDMS-DCDMS are shown in Fig.O Again it can be 
seen that the stick theory overestimates the magnitude 
of the force at small separations on both extension and 
retraction. In contrast, the slip theory with 6 = 2nm 
fits the measured results on extension down to several 
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FIG. 5: Vertical cantilever deflection Zc versus separation h 
at a drive rate of ip = 50/ims"^ for DCDMS-DCDMS (see 
Table [I}. The effective cantilever spring constant is fccff = 
2.22 N/m. Tfie circles are the analysed experimental mea- 
surement (every tenth point plotted), the dotted curve is the 
stick calculation, h — Onm, and the dashed curve is the slip 
calculation, b — 2 nm. The source of the raw measured data 
is Ref. 3. 
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FIG. 6: Deflection Zc at first contact versus drive velocity 
for DCDMS-DCDMS. The symbols are the analysed mea- 
sured data (raw data from Ref. 0), the dotted curve is the 
stick calculation b — Onm and a Hamaker constant of A — 
5 X 10"^" Jm~^, the dashed curve is the slip calculation with 
b = 2 nm and A — 5x 10"^" J ra~^ , and the dash-dotted curve 
is the slip calculation with 6 = 2 nm and yl = 1 x lO"^'' J m~^. 
The other parameters are as in Table HI The bars are a crude 
estimate of the error arising from the uncertainty in identify- 
ing first contact. 



FIG. 7: Adhesion defiection Zc versus drive velocity for 
DCDMS-DCDMS. The symbols, curves, and parameters are 
as in the preceding figure. The adhesive tension is = keszc. 



give an idea of the sensitivity to the Hamaker constant. 
Most of the measured data lie between b = Onm and 
6 = 2 nm. It is a little unexpected that the slip length for 
the hydrophobic low energy surface DCDMS should be 
less than the slip length for the hydrophilic high energy 
surface Si. 

Figure [7] shows the adhesion deflection as a function 
of drive velocity, again including repeat measurements. 
Again most of the measured data lie between 6 = nm 
and 6 = 2 nm. The Hamaker constant that was used, 
A = 5 X 10~^° Jm~^, is forty times larger that was used 
for Si-Si. It can be seen that a Hamaker constant of 
yl = 1 X 10^-^" J m^^ significantly underestimates the ad- 
hesion. (One might argue from the low velocity data 
that a value of A = 4 x 10~^° J would give a bet- 
ter fit. One might also argue that a smaller slip length 
would better fit the high velocity data.) To be con- 
crete, at a drive velocity of ip = 50^ms~^, the mea- 
sured adhesion deflection is 134.7 nm. For 6 = 2nm and 
A — 5x 10~^° J m~^, the calculated adhesion is 133.8 nm, 
and for 6 = 2nm and ^ = 1 x 10~^° Jm~^, it is 81.3nm. 
These calculations were performed with the variable drag 
algorithm. Using constant drag instead, for 6 = 2 nm and 
A = 5x 10~^° J m~^, the calculated adhesion is 134.9 nm, 
which is a minor change. (In this case the deflection at 
first contact increases by 0.9 nm using constant rather 
than variable drag.) 



nanometers from first contact. It is also in surprisingly 
good agreement for the measured retraction force. 

The extension deflection at flrst contact is shown in 
Fig. [5] as a function of drive velocity, including a number 
of repeat measurements. There appears to be a larger 
scatter in the measured data in this case compared to Si- 
Si analysed in Fig. |31 The crude error estimate appears 
to be on the low side. Calculated results for A = 5 x 
10"^" Jm~^, and for A = 1 x 10~^° Jm~^ are shown to 



C. Si-OTS 

Figure [8] shows the measured and calculated drainage 
force for Si-OTS. Given the similarity in contact angles, 
the OTS surface is expected to exhibit similar slip char- 
acteristics to the DCDMS surface, 6dcdms = 2 nm found 
above. Although 6ots = 2nm is not a bad fit, it can be 
seen that at small separations on extension, at this veloc- 
ity the measured data is better fitted with a larger slip 
length, 6oTS = 10 nm. Both curves use 6gi = 3nm, which 
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FIG. 8: Vertical cantilever deflection Zc versus separation h 
at a drive rate of ip = 50/ims~^ for Si-OTS (see Table |l|. 
The effective cantilever spring constant is fcefi = 1.34 N/m. 
The circles are the analysed experimental measurement (every 
tenth point plotted), the dotted curve is the stick calculation, 
6si ~ boTS = Onm, the dashed curve is a slip calculation, 
bsi — 3 nm and 6ots ~ 2 nm, both with a Hamaker constant 
of j4 = 1 X 10^^" Jm~^, and the dash-dotted curve is a slip 
calculation with fesi = 3nm, 6ots = 10 nm, and A = 2 x 
10"^" Jm~^. The source of the raw measured data is Ref. 0. 
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FIG. 9: Deflection Zc at first contact versus drive velocity for 
Si-OTS. The symbols, curves, and parameters are as in the 
preceding figure. The bars are a crude estimate of the error 
arising from the uncertainty in identifying first contact. 



FIG. 10: Adhesion deflection Zc versus drive velocity for Si- 
OTS. The symbols, curves, and parameters are as in the pre- 
ceding figure. The adhesive tension is Fz = kesZc- 



As for the preceding cases, the adhesion increases with 
increasing velocity, and at each velocity it has about the 
same magnitude as the deflection at first contact. It can 
be seen from Fig. [10] that the calculation 6si = 3 nm, 
^OTS = 10 nm and A — 2 x 10~^°Jm~^ overestimates 
the adhesion by about 150% at low drive velocities and 
underestimates it by about 10% at high velocities. De- 
creasing the Hamaker constant to A = 1 x lO"^'^ J m~^ 
improves the fit at low velocities, where it overestimates 
the adhesion by about 40%, but makes it worse at high 
velocities, where it underestimates the adhesion by about 
15%. Overall, the slip length 5ots = 2nm is in signifi- 
cantly better agreement with the measured adhesion than 
is the slip length 5ots = 10 nm. 

For the case of Si-OTS, the defiection at first contact 
cannot definitively distinguish between feoTS = 2 nm and 
^OTS — 10 nm. However, based on the results for the 
adhesion defiection, one can say that the slip length is 
much closer to 6ots = 2nm than it is to 6ots = 10 nm. 
The former value is the same as that found for the similar 
hydrophobic surface, 6dcdms = 2 nm. 



III. CONCLUSION 



was the value found above. For the retract curve, both 
slip lengths are about equally good. The calculations in 
the figure were obtained with the variable drag algorithm. 
There is negligible difference for constant drag. 

For the deflection at first contact. Fig. ^ as expected 
the stick theory overestimates the repulsion due to the 
drainage force. Within the experimental scatter, there 
is little to choose between a slip length of 6ots — 2 nm 
and &OTS = 10 nm. Perhaps the former is slightly better 
at high velocities and the latter is slightly better at low 
velocities. The Hamaker constant has almost no effect 
on extension (approach). 

The adhesion defiection for Si-OTS is shown in Fig.lTUl 



The non-linear algorithm for the analysis of atomic 
force microscopy data has been improved and made more 
reliable here. Specifically, the conversion of voltage to 
cantilever deflection is now more reliable on extension out 
of contact. Also the spring constant is more reliable as it 
is now determined by a least squares fltting procedure. 

Previously measured raw data^^i^ was re-analyzed with 
the non-linear algorithm. For silica, the slip length was 
found to be 3 nm, for dichlorodimethylsilane it was 2 nm, 
and for octadecyltrichlorosilane it was 2 nm. These val- 
ues are from 3-15 times smaller than those originally ob- 
tained with the linear analysis of the same data.'2iS ft is 
of note that the two low energy surfaces have a smaller 
slip length than the high energy silica surface. 
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Hamaker constants for the van der Waals force were 
also obtained from the measured retraction data. For 
Si-DOPC-Si the Hamaker constant was found to be 
1 X 10-21 J for DCDMS-DOPC-DCDMS it was 

5 X 10-20 Jni-2, and for Si-DOPC-OTS it was 1 x 
10-20 J ni-2. In so far as OTS is similar to DCDMS, 
the fact that the asymmetric system lies intermediate 
between the two symmetric systems is physically reason- 
able. 



Acknowledgments 



The raw data used here and originally analyzed in 
Refs Hi were measured by Liwen Zhu under the super- 
vision of Chiara Neto, and I thank Dr Zhu for providing 
the data. 



^ Neto, C, Evans, D. R., Bonaccurso, E., Butt, H. J., and 
Craig, V. S. J. (2005), Rep. Prog. Phys. 68, 2859. 

^ Bocquet, L. and Charlaix, E. (2010), Chem. Soc. Rev. 39, 
1073. 

^ Zhu, L., Attard, P., and Neto, C. (2012), Langmuir, 28, 
7768. 

* Attard, P. (2012), larX"iv:1212.3019t r2 . 

^ Zhu, L., Attard, P., and Neto, C. (2011), Langmuir, 27, 
6701. 

^ Zhu, L., Attard, P., and Neto, C. (2011), Langmuir, 27, 
6712. 

Zhu, L., Neto, C. and Attard, P. (2012), Langmuir, 28, 
3465. 

* Craig, V. S. J. and Neto, C. (2001), Langmuir 17, 6018. 

^ Attard, P., Pettersson, T., and Rutland, M. W. (2006), 

Rev. Sci. lustrum. 77, 116110. 
^° Higgins, J., Proksch, R., Sader, J. E., Polcik, M., McEn- 

doo, S., Cleveland, J. P., and Jarvis, S. P. (2006), Rev. Sci. 

lustrum. 77, 013701 . 
" Attard, P., Schulz, J., and Rutland, M. W. (1998), Rev. 

Sci. lustrum. 69, 3852. 
12 Attard, P., Carambassis, A., and Rutland, M. W. (1999), 

Langmuir 15, 553. 
" Vinogradova, O. I. (1995), Langmuir, 11, 2213. 

Attard P. and Parker J. L. (1992), Phys. Rev. A 46, 7959. 

Erratum, (1994) Phys. Rev. E 50, 5145. 
1^ Southwell, R. V. (1936), An Introduction to the Theory of 

Elasticity, (Oxford, London). 
1^ Vinogradova, O. L and Yakubov, G. E. (2003), Langmuir, 

19, 1227. 



Appendix A: Improvements in Atomic Force 
Microscopy 

One of the crucial issues is the calibration of the can- 
tilever that is used to measure the surface forces: us- 
ing too small a value for the spring constant causes the 
slip length to be overestimated and conversely if it is too 
large. Hence motivated by the need for accurate values, 
one significant improvement in atomic force measurement 
methodology was the in situ measurement of the spring 
constant of the cantilever.— This was done by measuring 
the drainage force at large separations, and only required 
values of the radius of the colloid probe and the viscos- 
ity of the liquid, both of which can be readily measured. 
This method has three advantages: First, as an in situ 
measurement the effective spring constant that emerges 



is the one that is required for the actual measurements 
that are being performed. In particular, it automatically 
takes into account the position at which the colloid probe 
is glued to the cantilever (the cantilever spring constant 
varies as the cube of the distance from the base), and also 
the tilt of the cantilever, which is typically about —11° 
(the force measuring spring constant varies as the square 
of the cosine of the tilt angle). Second, the statistical 
error is much reduced, being on the order of 1% for the 
hydrodynamic method compared to 10-20% for the ther- 
mal method, for example. Third, as has been pointed 
outfit the most common conventional method of spring 
constant calibration, the thermal method, suffers from a 
mathematical error in an early publication,^" which has 
been faithfully reproduced in the built-in software of at 
least one brand of atomic force microscope, and which 
causes a systematic overestimate of the spring constant 
of 15%-30%<2ii^ (The correct formulae for the thermal 
calibration is given in Eq. (35) of Ref.i.) 

It is again emphasized that an error of 10% in a first 
order quantity such as the spring constant can translate 
into an order of magnitude error in a second order quan- 
tity such as the slip length. 

If the spring constant calibration can be taken as the 
first generation, then the second generation of improve- 
ments to slip length and drainage force measurement 
were undertaken by the present author with colleagues 
Zhu and Neto and presented in a series of papers.^^"— 
Briefly these are: 

• a better numerical method of calculating the 
drainage force that is independent of the measure- 
ment, and a way of fitting it to the measured data 
that is more sensitive to the slip length than previ- 
ous methods 

• in situ temperature measurement, which is impor- 
tant because the viscosity is sensitive to tempera- 
ture, which can vary over the course of a series of 
measurements 

• a 'blind test' protocol, which is used to establish 
the statistical accuracy of the methods used to fit 
the spring constant and the slip length 

• care to exclude particle contamination, which cre- 
ates an artefact that shows up as a very large slip 
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length 

• an account of thermal drift and virtual deflection in 
the data analysis, which can effect both the spring 
constant determination and the slip length fit 

• accounting for the variation in the drag force on 
the cantilever with cantilever deflection, the ne- 
glect of which causes the slip length to be over- 
estimated and to become cantilever- and spring 
constant-dependent 

Of these various improvements, perhaps the most novel 
is the variation in the drag force with deflection. This 
explained quantitatively many of the puzzling results for 
the slip length that had been reported previously in the 
literature, including the observation that it appeared to 
be larger for softer cantilevers and that it appeared to 
depend upon the type of cantilever used. Whether or 
not the non-constant cantilever drag force is important 
in other measurements with the atomic force microscope 
depends upon experimental details. In general terms, 
reducing the drive speed and increasing the stiffness of 
the cantilever reduces the influence of the variation in 
drag force. 

What might be called the third generation of improve- 
ments consist primarily of taking into account non-linear 
effects in the raw atomic force microscope data4 Such 
non-linearity is a general feature of the atomic force mi- 
croscope, and the corrections for it have general applica- 
tion beyond the drainage force. The non-linearity is illus- 
trated in Fig. [U where the measured photo-diode voltage 
is plotted against the piezo-drive position approaching 
and in contact. The conventional analysis of atomic mi- 
croscope force data begins by taking the slope of the line 
in contact, V/nm, the so-called constant compliance fac- 
tor, and using it directly to convert a change in voltage 
to a change in cantilever deflection, and hence, using the 
spring constant, to a change in force. This is how forces 
are measured with the atomic force microscope. However 
as can be seen in Fig. [TJ the contact region does not have 
constant slope, and whatever value of the slope is chosen 
as the calibration factor has a quantitative effect on the 
data analysis that gives the measured forces. The effect 
is quite dramatic in the contact region, as can be seen in 
the inset. However, even though the effect seems small 
for the pre-contact drainage force on the scale of the flg- 
ure, it can change by more than an order of magnitude 
the fltted slip length^ 

In Ref. y two sources of non-linearity were anal- 
ysed: non-linear photo-diode response and non-linear 
cantilever deflection. It was concluded that in typical 
experiments the latter was negligible, and that it was 
the photo-diode that was responsible for the non-linearity 
evident in Fig. [TJ (Most likely, it is the elliptical cross- 
section and Gaussean intensity distribution of the light 
beam that directly causes the non-linearity.) The spread- 
sheet algorithm that was developed in Ref. to analyze 
surface forces with a non-linear contact region has a num- 
ber of features: 



• it removes the ambiguity in the calibration factor 
by using a polynomial fit for the non-linear conver- 
sion of the measured voltage to cantilever deflection 



• relying upon earlier analysis jiiii^ it takes into ac- 
count cantilever tilt, friction, and torque to distin- 
guish between the intrinsic and the effective can- 
tilever spring constant, and to obtain the actual 
pre-contact force from the raw voltage 

• it automates the identiflcation of the zero of sep- 
aration, which avoids human intervention or bias, 
and which can be particularly problematic when 
the linear analysis gives an unphysical contact re- 
gion 



• it allows the simultaneous analysis of both the ex- 
tension and the retraction data so that the drainage 
adhesion can be used as an additional constraint on 
the determination of the slip length 

This third generation analysis was applied to the same 
raw data as was previously analysed with the second 
generation technique, namely the drainage force between 
silica surfaces in di-n-octylphthalate (DOPC). Whereas 
the second generation linear analysis had given the slip 
length as 11±2.5 nm, - with the non-linear analysis it was 
found to be 3 ± 1 nm.** In addition, the linear analysis 
had shown that the slip length decreased with increasing 
shear rate, whereas the non-linear analysis had shown 
that it was constant. 

The reasons why the linear analysis had these partic- 
ular quantitative and qualitative failings was discussed 
in detail in the conclusion of Ref. 0. Briefly, due to 
the fact that the compliance curve. Fig. [1] is concave 
down, the compliance slope taken from the tangent at 
first contact is an underestimate of the factor required 
in the non-contact region. This means that the effective 
spring constant fitted at large separations is too small 
{kcs — 1.5 N/mf^ compared to fceff = 1.68 N/mi obtained 
with the non-linear analysis), and consequently the slip 
length required to fit the data at intermediate separations 
is too large. This overcorrects the problem closer to con- 
tact, where the calibration factor that is used is closer to 
the real one. Since the shear rate increases with decreas- 
ing separation, the reduced slip length required at small 
separations was interpreted as being due to shear rate 
dependence.^ In the non-linear analysis of Ref. 4, a sin- 
gle slip length gave a good fit over all separations (down 
to about 1 nm) and for all drive velocities. As well, the 
retraction data, which was not originally analysed^, was 
also well fitted using the non-linear analysis, over almost 
the whole separation regime and at all drive velocities 
using the same slip length. 
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Appendix B: Non-Linear Data Analysis Algorithm 
for a Spread-Sheet 



This section sets out in order the steps for the spread- 
sheet algorithm for the non-hnear analysis of raw data 
measured with the atomic force microscope. The section 
is written in the form of a recipe, and only limited jus- 
tification and explanation is given. The formula, with 
small modifications, are derived in Ref. 0, with the linear 
cantilever coefficients originally derived in Refs. [TTIIT2I . 

Here it is assumed that the raw extension and retrac- 
tion data are separated. These quantities are denoted 
'ext' and 'ret' respectively, and each is analysed inde- 
pendently except as noted below. Further 'c' refers to 
contact, 'nc' refers to non-contact, and 'b' refers to base 
line. 

In the atomic force microscope modeled here, the 
piezo-drive is connected to the cantilever holder and is 
above the substrate. In extension the velocity of the 
piezo-drive is negative, Zp < 0, and in contact the move- 
ment of the contact point is equal and opposite to that 
of the piezo-drive, Azc = — Azp. In other models of the 
atomic force microscope, it is the substrate that is con- 
nected to the piezo-drive, so that ip > on extension and 
in contact Azc = Azp. These other models may be anal- 
ysed with the formulae given in this paper by negating 
the piezo-drive position in the raw data, Zp —Zp. 



1. Base Line 

Select the base line region, where the raw voltage is a 
linear function of the drive distance. Do not include the 
initial extension data or the final retraction data where 
the piezo-drive is not moving at constant velocity. Do 
not include (or include as little as possible) of the region 
where the data is clearly curved due to the drainage force. 

Independently fit two straight lines to the measured 
voltages in the base line regions: 



y'GXt 
b 



(^p) 



and 



cxtl 



'pbJ 



(Bl) 



(B2) 



2. Force Asymptote 

Calculate the linear asymptote to the drainage force in 
the base line region: 



Fr(^p)=Fr 



f; 



and 



This requires that the separation be known, h 



, cxt/rct 

+ ^0 • 



The constants 



-0 



and 



-0 



(B3) 



(B4) 



Zp + 



are obtained in 



"^pb J; 



rot] 
^pbJ- 



TO 61 below. The vertical deflection of the contact point 
is Zc — Fz/kcs- With C = —GnriR^, one can set up the 
iteration 



F. 



(0) 



Cz, 



zq + Zp 



and Fi") 



Cz 



zo ■ 



(n-1) ■ 



(B5) 

Generally only zero or one iterations are required. One 
can recognize here ft.^"-* = zq + Zp + k~f^Fz"\ For the 
base line, choose Zp = Zph and this iterative formula gives 
the base line constants Fb and = —F^/h. Obviously 
one uses either 2°^^' or ip°* for extension or retraction, 
respectively. 

Its worth mentioning that the algorithm given in i)C 71 
below was also used for the extension branch. It is 
slightly more robust than the present one, but it makes 
no difference to the fitted spring constant. The present 
algorithm has the advantage that it does not need to cal- 
culate the retraction curve from contact. 

In the above, 77 is the viscosity, R is the radius of the 
colloid probe, and kos is the effective cantilever spring 
constant. The method of determining the spring con- 
stant is given below. Fortunately, the results are not 
very sensitive to this value, and so in the first instance 
any estimate can be used, and this can later be replaced 
by the accurate value. 



3. Contact Fit 

Perform non-linear fits to the measured voltage in the 
contact regions: 

z^fiV) = 0'=''* + W^'^V + c'^'^V^ + d^'^V^ + . . . , (B6) 

and 



These include the effects of virtual deflection, thermal 
drift, drag, and drainage force at long range. Different 
spread sheet programs fit straight lines in different for- 
mats. It is assumed that the user can convert the fit- 
ted form to the above format, with Zph chosen to be in 
the middle of the selected base line region by the user. 
Experience indicates that best results are obtained with 



■^pb 



"pb- 



PC 



{V) 



(B7) 



Note the distinction between Zpc, which is the piezo-drive 
position at a given voltage when the probe is in contact 
with the substrate, and Zc, which is the vertical deflection 
of the contact position (apex) of the probe. In practice 
four terms were found adequate. One should probably 
avoid increasing the number of terms. 



10 



It is important to appreciate that the extension fit only 
covers voltages strictly greater than the voltages mea- 
sured in the non-contact extension region. In contrast, 
the retraction fit covers almost all measured voltages, 
both extension and retraction, contact and non-contact. 
(The exception is the voltages around the minimum of 
the retract curve, which is not a large extrapolation any- 
way.) It has been found, for example, that the slope of 
the contact voltage evaluated at the extension base line 
voltage is unreliable when extrapolated from the exten- 
sion fit, but is reliable when interpolated from the re- 
traction fit. (Reliability here was judged by whether or 
not it remained more or less constant over a long series 
of measurements.) The difference was around 5% in a 
typical case. For this reason the retraction fit will pri- 
marily be used for the non-contact data in what follows. 
This is one of two improvements on the algorithm given 
in Ref. 4, 

It is not appropriate to use the retraction fit for the ex- 
tension data in contact. This is because friction plays an 
opposite role in the two cases, and also because in contact 
the extension fit is strictly an interpolation for extension 
in contact (except possibly near the ends). Hence as ex- 
plained below a switch is used to apply the two different 
fits to the extension data in the two regimes. (The re- 
traction fit can be used for all the retraction data.) Such 
a switch is only necessary if one wants to look at the 
topography in contact. It does require a fairly accurate 
determination of the point of first contact. 



4. Compliance Slope 



The retraction gradient is 



dZn 



(B9) 



These slopes are dominated by the change in angle of 
the cantilever due to the surface forces in contact. How- 
ever, they also contain a contribution from the base line 
slope, which should really be removed. One could define 



'CXt ocxt 

cbf = Pcb 



T/CXt 



OGXtjL-i pox 

Pcb '^cS 



(BIO) 



and similarly for f3"'^ 



^^f. (Note that in contact, dzc = 
— dzp.) These corrections are generally negligible. In 
practice this correction was implemented with = 
(to avoid circularity) even though the change was typi- 
cally only 0.05%. 

For completeness, the rate of change of voltage with 
vertical contact position out of contact is 



dV dV 


d9 








d^c^ ~d9 
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dzc 


d9 


dZp 


ret 
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ret d^^e 



^ret 2c™ty -f ZdF^^V^ 



(Bll) 



This holds for both extension and retraction out of con- 
tact. 



One eventually has to convert the base line voltages 
to a deflection angle, and for this the compliance factor 
derived from the gradient of the contact fit evaluated at 
the base line voltage is required. The constant compli- 
ance gradients for the base line are defined as 



5. Contact Position 

The initial estimate of the contact position is where 
the contact voltage equals the constant part of the base 
line voltage. 



/□ext 
Pcb 



dVr\zp) 



dV 




d^ 
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T/"ext 
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oxt d^:p 


c d^c 


ret 
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dVr'izp) 



Q,cxt/^rct 



(B8) 



Here 9 is the angle of defiection of the cantilever. The 
formula for a = d9/dzc is given in TO 91 

Notice that only the retraction fit is used for this. In 
order to convert from the extension slope to the retrac- 
tion slope one has to account for difference between the 
rates of change of angle on extension and retraction in 
contact, which is the reason for the ratio in the numera- 
tor. 



and 



ext 
^pcb 



'pcb 



^retyjoxt , 



(B12) 



(B13) 



Again the retraction fit is used for both, but they are 
evaluated at the respective base line voltages. 

Now this is corrected so that the fitted voltage equals 
the actual linear base line voltage at that position, 
Vb* = ^pfiV^'^'i^T))- A Taylor expansion yields 

evt, * ^vf dZr 



'pcb 



^pcb 



( /OCXt \ — 1 T/CXt' 

IPcb ) 



TACXt^.^CXt,*\ 

^b (^pcb ) 



-pcb 



V'CXt 
b 



'pb 



l^ext ^ext 
Pcb ^pcb 



^pb 



'■pcb 



/Qext 
Pcb 



T/ext' 



,(B14) 



(B15) 
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One has an analogous result for retraction, 



ret,* ^cb pcb 

■^pcb ~ Qvet 
Pch ' 



T/ret' ret 
•^b -'-pb 



(B16) 



In so far as the compliance slope is generally much 
greater than the base line slope, the difference between 
^pcb* -^pcb generally small, 1.5 nm in a typical case. 
This is negligible at large separations but is important at 
small, particularly if one wants the topography in con- 
tact. 



6. Zero of Separation 

The deflection of the contact position due to the force 
in the base line region is for extension, 

^cb*(%) = ^efl ■'^b''*(^p) 

= {f^""' + Fb""'' [zp - z;f] } , (B17) 
and for retraction, 



zlbizp) — k^g F^"^*" (zp) 



^pbj 



} ■ (B18) 



These linear results are defined as holding until con- 
tact, even thoiigh the actiial measured forces depart sig- 
nificantly from the base; line. Recall that the contact 
point is defined as the intersection of the extrapolated 
base line and the extrapolated contact curve. 

The separation is in general h = Zp + Zc + zq, with the 
constant zq to be now determined. At contact, Zp = z*^^^, 
one must have h = 0. Hence 



~ext _ cxt,* 
^0 - ~^pcb 



cxt/ cxt,* N 

^cb (Zpch ) 



and at z. 



pcb 



ret _ _ ret,* 



'0 



= -z. 



pcb 



_ ^ret/ ret,*x 
^cb l^pcb ) 



(B19) 



(B20) 



These give the zero of separation. 

It will be noted that there is no direct human inter- 
vention in the choice of the zero of separation. These 
constants emerge automatically from the flts to the base 
line and contact voltages. 



7. Base Line Deflection Angle 



The base line angle for extension is 



dzc 
dV 



Zp 2pb J 



a 



next 

,,ext 



Aext 



/3ext 
Pch 



V^-''[zp-z;f]. (B21) 



The two final equalities follow either by direct substi- 
tution, or else by evaluating the derivatives in the first 
equality in contact on extension. Similarly, that for re- 
traction is 



qrct 
■'b 



^ T/ret' r _ 

oret "^b L^P ^1 
'^cb 



ret] 
pbj 



(B22) 



There are two undetermined constant angles here: 
and 9l^*'. Since the difference in angles is linearly propor- 
tional to the difference in base line voltages, a symmetric 
choice is 

" [V;-* - V^^'] , (B23) 



oret /lext 

— ~''h — 



with Vh = [V^"^ + V^''^]/2. In fact, these two constants 
drop out of the final formula below and so the actual 
values are immaterial. 

The deflection angle due to the base line force for re- 
traction is 

^"f*^^'') = £ {^"* + ~ } ' ^^^^^ 

and that for extension is 



9 cxt 
bf 



^^p) = ^ {^^^' + " 4^'] } ■ ^^^^^ 



For convenience below the constant parts of these may 

'bf — "^b 



be defined as Olf = aF{;^/kes and 6'°f = aF^^^yk^s. 



8. Deflection Angle 

All the formulae above have been for constants, or lin- 
ear base line fits, or polynomial contact fits. Now is given 
the formulae that convert the measured raw voltage into 
a deflection angle of the cantilever. Obviously the formu- 
lae will invoke the fits and constants given above. 

In the atomic force microscope one has a set of data 
pairs {zp, V}, separately for extension and for retraction. 
The aim is to convert each data pair to separation and 
force. Since the latter two are linear functions of the de- 
flection angle, and it is the components of the angle that 
arc linearly additive (even though the voltage is a non- 
linear function of the angle), the formula are obtained for 
the deflection angle. 

The deflection angle on the retract contact curve is 



a 



'pcb 



(B26) 



The constant 61^*^'* is the deflection angle at the inter- 
section of the contact curve and the linear base line, and 
must be the same for both (since by definition they have 
the same voltage at this point). 
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The total angle on the retract contact curve is 

Clc(^) = + env) + C(4c'(^)) - €fKciv))- 

(B28) 

The constant tilt angle is ^o- (All angles are in radians; 
the tilt angle is negative.) Everything else on the right 
hand side was defined above. 

The total angle at a given voltage and a given position 
on extension is the sum of the tilt angle, the deflection 
angle, and the base line angle less the base line deflection 
angle, 

etTdv, zp) ^Oo + e^-Hv, z,) + (^p) - (^pX-B29) 

Since there is a one to one correspondence between the 
photodiode voltage and the total cantilever angle, this 
must equal the total angle in contact on retraction at 
the same voltage, 9'^^^{V,Zp) = 0[ot.c(^)- Hence the de- 
flection angle at a given voltage and a given position on 
extension is 

e^-\v, zp) = Cc(^) - ^0 - oi-\zp) + eifizp) 
= env) + - eif{z;i\v)) 

-ei-\zp) + e^^fiz,,). (B30) 

Since 9'^^^{V^^^ , z^^) = 9^^, which was give above, the 
constant terms may be collected and this expression may 
be rearranged as 



for non-contact extension, because only interpolation is 
required. In the non-contact regime, friction has no in- 
fluence, but in contact it does, and this approach messes 
up the angle for extension in contact. In the usual case 
there is no need for data in contact. However for com- 
pleteness, the formula to be used for extension in contact 
is as above with 'ret' change to 'ext'. 



'iV,Zp)^ 

-a-* [z;f{v) - z;%\ - — ^T/r Kfi^) 

Pch 

- ^J^b"' Hfiv) ~ z;^i] + ^vr' [z, 
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(B33) 



ext] 
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^^^b L^P ^pb J + 
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bf ■ 



9. Cantilever Characteristics 

For a cantilever of length Lq (more precisely, the probe 
is attached a distance Lq from the base of the cantilever), 
inclined at an angle to the horizontal of 6*0 < 0, with col- 
loid probe of radius R and diameter L2, (more precisely, 
L2 = R\/2~rYcos6o) , one can define Co = cos^o and 
Sq = sin 610. The ratio of change in deflection angle to 
the change in vertical deflection in the linear cantilever 
regime isiiiiii^ 



-a" 



(B31) 

xt ] 
pcbj 



^F-^' [z-(y) - z-*] + ^vr' [z, - ^rb' 



a 
a 



^pb 



^bf ; 



smce 



^pc (^b )• This result hold on extension out 
of contact. (See the end of this section for the contact 
formula.) 

The analogous result for retraction is 



'iV,Zp) 



(B32) 
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fccff 

a 

fccff 



2pb 



bf 



From these deflection angles, the vertical deflection of 
the contact position can be obtained, Zc = a^^0, with 
the constant given in fjBjO] From the deflection and the 
constants given above, the separation can be obtained, 

h ^ Zp + Zc + Zq. 

As mentioned above, it is most accurate to use the re- 
traction contact voltage flt to obtain the deflection angle 



D{ii)Co + L2S0 



R Sq9o 



+ L2Coeo} (B34) 



The friction coefficient is positive, ^ > 0. There are 
three values of the ratio: extension in contact, retraction 
in contact, and non-contact. These are denoted 

EE a(^), a"' = a(-^), and a = a(0). (B35) 

The ratio of the deflection to the deflection angle that 
appears here is 



(B36) 



10. Spring Constant Determination 

As discussed previously;^ the effective spring constant 
fceff and the intrinsic spring constant are related by 



^eff 



2Lo/3 



Co -I- 2L2Sq/Lq 



fco 



(B37) 
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The effective spring constant results after taking into ac- 
count the cantilever tilt and the torque on the colloid 
probe. It gives directly the force acting on the probe 
from the vertical deflection of the contact point, 



Fz = keSZc 



(B38) 



(The vertical deflection of the contact point arises in the 

atomic force microscope force measurement in the con- 
version factor for voltage, because the calibration of this 
is based on the fact that the vertical deflection of the 
contact point is equal and opposite to the movemc^nt of 
the piezo-drive when the probe is in contact with the 
substrate.) The intrinsic spring constant is a material 
property of the cantilever and is what is measured in, for 
example, the thermal calibration method. 

The results given above convert the raw atomic force 
microscope voltage into vertical cantilever deflection, so 
that one has the data pairs {zp,Zc} for both extension 
and retraction. The value of the spring constant only 
entered this conversion in the calculation of the base line 
force Fh{zp). This has only a weak influence on the re- 
sults. So what one does in practice is first estimate an 
initial value /c^g , then use the following procedure to fit 
the effective spring constant to the {^p, Zc} data. One can 
use this fitted value to replace the initial value used for 
the base line force, and then repeat the fit. In practice, 
repeating the fit made negligible change. 

To fit the effective spring constant the iterative expres- 
sion for the force given above is used but for the deflection 
rather than the force, 



^(0) ^ k^Czp^ 

ZO + Zp 



Zo + Zp + Zc 



(n-1) 



(B39) 

In practice two iterates were used, which means three 
columns of data each for extension and retraction, with 
each row corresponding to a given measured Zp. Note 
that this is the stick drainage force formula; the slip 
length does not come into it. 

It should be mentioned that this ignores the rate of 
change of the cantilever deflection, Zc- This is valid at 
large separations where \zc\ <C |ip|. 

The calculated deflection was converted to voltage by 
using 

y-'(zp) = [KsF^'^'i^p) - ^c] + Vnzp), (B40) 

and similarly for retraction. The linear base line force 
that appears here is the one with the original estimate 
of the effective spring constant. The use of the base line 
calibration factor /3cb means that this is restricted to the 
linear regime, which means that the deflections cannot 
be too large, which again restricts the formula to large 
separations. The square of the difference between this 
voltage and the measured raw voltage was summed over 
the specified ranges for extension and retraction and min- 
imised to obtain fceff- 



The iteration procedure fails for large deflections. De- 
pending upon the system, large can mean 1 10 nm. It is 
clear when it fails, because the iterations do not converge. 
Care has to be taken to exclude failed Zc from the error 
estimate. The way this was done is to specify a range 
|-_2,smaii^ ^largej Q^gj. -^^jjjQjj |;}^g errors were summed. The 

upper limit is close to the start where the piezo-drive is 
moving at uniform velocity; it is typically the same at the 
base line upper limit. The lower limit is greater than the 
first position when the iteration procedure fails. Typi- 
cally a range of 3-5 /zm was used for the estimate of the 
error in the fit. One can always tell if the range is appro- 
priate and if the fit is good from a plot of the measured 
and the calculated deflections. 



Appendix C: Variable Drag Derivation and 
Algorithm 

1. Drag Length 

At large separations the drag force on the cantilever is 
constant and can be described by an effective drag length. 



pext 
drag 



-eTTTyi^'^'Ldrag 



(Cl) 



The effective drag length is typically a fraction of the 
length of the cantilever, but this obviously varies between 
cantilevers. 

The constant drag force is removed from the measured 

experimental data by the treatment of the base line dis- 
cussed above. Hence in most cases one does not need to 
know directly the drag force or the effective drag length. 
However in some cases, particiilarly for soft cantilevers, 
highly viscous liquids, and high driving velocities, the 
variation in the cantilever drag force with cantilever de- 
flection cannot be neglected. In this case the effective 
drag length is needed in order to calculate the variable 
drag force, as is discussed shortly. 

The difference in the base line voltage on extend and 
retract is due in part to the change in sign of the con- 
stant drag force, in part due to the change in sign of the 
drainage force in the base line region, and in part due 
to thermal drift. The latter two have to be subtracted 
in order to obtain the drag force alone. The change in 
voltage due only to thermal drift is 



ye^t' _ yret' 
Tiext' T?ret' 



-[' 

X [22^p,turn ~ •Zpb ~ •^pb] • 



(02) 



where ^;p,turn is the furthest position in contact, where the 
piezo-drive turns from extension to retraction, and Vb = 
^yext _|_ y^ctj^2. Clearly the contribution to the change in 
voltage from the drainage force is being subtracted here. 

This expression assumes that the piezo-drive velocity is 
constant on each branch. Since in fact it slows just before 
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the turn point, and it accelerates just after, it would be 
better to write this in terms of the change in time (since 
the thermal drift is constant in time). In practice this 
was found to make negligible difference in the drag force. 

With the above expression for the change in voltage 
due to drift, the drag force is 



pcxt 
drag 



^cff 



therm 



2A(Fb) 



(C3) 



R 



R 




where again the drainage force is subtracted. In prac- 
tice the term AVthcrm is often negligible. From this the 
effective drag length can be obtained. The standard de- 
viation in the effective drag length over all velocities is 
about 3%. 

There is an alternative way of obtaining the effective 
drag length that also works well and serves as a useful 
check or replacement for the above procedure. At the be- 
ginning of the extension run, the piezo-drive decelerates 
from Zp'^' = to z'^'^ = — |ip| over typically a distance of 
several hundred nanometers. Similarly, at the end of the 
retraction run it accelerates from ip*'' = |ip| to ip°' = 0. 
Depending on the model of the atomic force microscope 
and the settings, this data is recorded and appears as 
a distinct region at the terminus of the extended flat 
plateaux that represent the drag and drainage forces for 
constant velocity. The change in voltage (equivalently 
force) from the beginning on extension, or end on retrac- 
tion to the plateaux values gives the drag force (after 
the asymptotic value of the drainage force has been sub- 
tracted). This second method is almost identical to the 
first method neglecting AVthcrm- 



2. Drainage and van der Waals Force 

The drainage force is 



Fdrmn(h,h) = v(h). 



(C4) 



For the case of stick, v{h) — 1. In this case the force 
is attributed to Taylor. For slip and symmetric surfaces 
each with slip length b on each surface, the Vinogradova 
factor is^^ 



v{h) 







In 






-( 








36 \ 











(C5) 



For the asymmetric case, with slip lengths hi and 62, one 
defines the relative asymmetry as fc = [62 — 61] /61 and 
also A = bi-\- 62, 



C± EE 26i 



2 + k±\/l + k + k'^ 



(C6) 



With these the factor is^ 

-2Ah 2h 

v{h) = 



FIG. 11: Cantilever geometry in the atomic force microscope 
(not to scale). 
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.(C7) 



In both case slip is negligible at large separations, v{h) — > 
as /i — >■ 00. 

The van der Waals force is taken to be^^ 



-AR 



4 I h 



(C8) 



Here A is the Hamaker constant. This includes a short 
ranged repulsive force, which is derived from a Lennard- 
Jones 6-12 potential. The location of the minimum in 
the potential is here taken to be — 0.53 nm. The van 
der Waals force in general has no effect on the drainage 
force or the slip length on approach (extension). It does 
however effect the adhesion, as is shown in the results 
in the text. For convenience the sum of the point forces 
acting on the probe will be denoted = -Fvdw + ^drain- 
For use below, the derivative of the van der Waals force 



dh ~ 3h^ 



and those for the drainage force are 



dF, 



drain 



dh 



v{h) _ v'jh) 



(C9) 



(CIO) 



where the prime denotes the derivative with respect to 
separation, and 



OF. 



dh 



v{h). 



(Cll) 



C+C- 



3. Cantilever Analysis 

This analysis follows closely earlier analysis of the 
bending of an atomic force microscope cantilever by At- 
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tard and co-workers ^i^^iiii^ In turn, that analysis is 
based on the classical equations of continuum elasticity^iS 

Consider a rectangular cantilever of length Lq and 
width w ^ Lq, tilted at an angle of 6*0 < to the hori- 
zontal with long axis in the j/z-plane. Define Cq = cos 
and 5*0 = sin^o- The tilt angle is small, but no lin- 
earisation will be performed. Angles are measured in a 
clockwise direction, and the y-coordinate is measured left 
to right. There is a spherical colloid probe of radius R 
attached to the end of the cantilever. There is a planar 
substrate in the a;y-plane located beneath the cantilever 
such that zero separation for the undeflected cantilever 
occurs when the base of the cantilever is at Zp = 0. 

Let (T e [0, Lq] denote the running coordinate along the 
cantilever, with cr = being the base and a = Lq being 
the end. Note that this runs right to left, the opposite 
sense to the y-coordinate. Let (^{a) denote the deflection 
orthogonal to the original axis, with = C(^) being 
the deflection of the end. The deflection angle is 9 {a) — 
dC((7)/d(T, or, equivalently. 



C(a) = / da' 



(C12) 



The boundary condition is C(0) = 9{0) = 0. The deflec- 
tion angle of the end, 9^ = 9{L), is what is measured 
in the atomic force microscope. The deflection will be 
assumed small and everything will be linearised with re- 
spect to it. 

In the laboratory frame, the undeflected cantilever is 
described by 



yo(o-) = -Cqcf, and zo(cr) = Zp + Saa. 



(C13) 



Here Zp is the piezo-drive position, which is the location 
of the base of the cantilever. 

The deflected cantilever has laboratory coordinates 

y{a) = - [ do-' cos(6'o +6'(cr')) 



and 



= -Co(7 + So / da' 9ia') + 0{9^) 
- yo{<j) + 5oC(a), (CM) 

zp+ f da' sm{9o+9{a')) 
Jo 

= zp + Soa + Co f da' 9{a') + 0(9^) 
Jo 

= zo(a)+CoC(a). (C15) 

With total angle ^tot ^ 9o + 9, a colloid probe of radius 
R has lever arm 



L 



2(,t7tot; 



= i?V2 + 2cos6itot 

= Ry/2 + 2Co - 2Sq9 + 0{9^) 



So9 



2[1 + Co] 
L,-^pl+0{9% 



0(9^) 



where L2 = -^2(^0) = RV^ + SCq. The change in the 
length of the lever arm with deflection angle is small. 

The change in the point of closest approach of the col- 
loid probe (contact position) is horizontally 

Vc = SoCl - L2{9tot) cos(6'tot)6'tot + ^2(^*0) cos(6'o)6'o 



and vertically 



L2C0 — L2So9q — Co9o 



R^Sq 



9l, (C17) 



CoCl + sin(0tot)i2(^tot)^tot — SqL29o 

R^Sq 



= CoCl + 



Sq + Co9q — Sq9o- 



^2 



L29l 



— CoCl + SqL29l. 



(C18) 



These are given to linear order in the deflection and rep- 
resent the change from the position of the undeflected 
cantilever. Note that to leading order the change in hor- 
izontal position has the opposite sign to the deflection 
and the deflection angle. This determines the sign of the 
friction force below. 
The separation is 



Zp -\- Zq 



(C19) 



As mentioned above, the base of the cantilever has been 
positioned relative to the substrate so that the unde- 
flected cantilever is in contact when Zp = 0. 

Let Fz and Fy be the components of the force act- 
ing acting on the colloid probe at the point of closest 
approach. The van der Waals and drainage force on the 
colloid probe are represented by F^, and the friction force, 
when present, is represented by Fy. Typically, Fy = ^F^ 
in contact on extension, where /i is the friction coefficient 
(see the above comment regarding the sign of the change 
in horizontal position) . These forces create a force on the 
end of the cantilever, 



Fl — FzC'o + FySo, 



(C20) 



which is orthogonal to the axis of the cantilever, and a 
turning moment (torque). 



Ml — L2FzSo — L2FyCo, 



(C21) 



which is positive in the clockwise direction. In these 
^'tot — 9o + 9 has been replaced by because only terms 
linear in the deflection are retained. 

Now add a force per unit length /(cr) in the can- 
tilever frame (i.e. acting orthogonal to the undeflected 
cantilever) . This is the drag (and possibly drainage) force 
on the cantilever, which will be given explicitly in the 
following section. In the cantilever frame, the turning 
moment about a due to the end force, end moment, and 
force per unit length is 



(C16) M{a)^[Lo-a]FL + ML- 



da' f{a')[a' ~a\. (C22) 
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From this the deflection of the cantilever isi^ 

C(ct) = r da" M{a")[a-a"] 

Jo 



Fl 
B 



2 

1 

B 



2B 







6 






. (C23) 
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Here B — EI is the beam elasticity parameter, which is 
related to the intrinsic spring constant of the cantilever 
by fco = W/Ll 

For the case of a uniform distributed force, f{a) — /, 



FlLI , LIMl , / 



B 3 2B 

Fl ^ 3Ml 

ko 2Loko 8fco 



B 



^0 



6 24 
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For the case that there is no drag or drainage force on 
the cantilever, /(cr) =0, the end deflection is 



Fl 
ko 



3Ml 



(C25) 



If in addition there is no lever arm, L2 = Ml = 0, nor 
friction force, Fy = 0, so that the end force is Fl — F^Cq, 
then the end deflection reduces to 



. FzCq FzCq 

Cl = — ; — , or Zc = — ; 

ko ko 



(C26) 



Here ClCq is the vertical deflection of the end of the 
cantilever in the laboratory frame. In this case with no 
lever arm, L2 = 0, the vertical deflection of the end is 
the same as the change in vertical position of the contact 
point, Zc = ClCq. 

This result does not agree with the result given by 
Vinogradova and Yakubov^^ for the same case (their 
Eq. (6), in the present notation, is Zc = F^ [3Co — l]/2fco) 
In the opinion of the present author, the error in Ref. [la 
arises in the very first equation where the authors mix up 
the internal coordinates of the cantilever with the lab- 
oratory coordinates. Consequently, the present author 
believes that all of the results in Ref. fie! are suspect. 

From the general expression for the deflection of the 
cantilever, the angular deflection of the end, Ol = 

dC(cr)/dcr|^^. , is 



2B 



LoMl 
B 



1 

B.10 
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This is important because it is what is actually measured 
in the atomic force microscope. For the case of a uniform 
force per unit length, /(ct) = /, this reduces to 



LIFl 
2B 



LqMl 
B 



6B ' 



(C28) 



In the case of only probe forces, (i.e. no distributed 
drag or drainage forces on the cantilever), a number of 
useful constants can be given. First define for conve- 
nience 



C+ = Co + fj.So and S- — Sq ~ /iCo, 



(C29) 



which differ from Cq and So only in contact since the 
effective friction coefficient is defined as 



/i, contact, extension, 
— /i, contact, retraction, 
0, non-contact. 



(C30) 



With these, the ratio of the end deflection to the surface 
force is 



F, 



9± 

ko 



3L2S- 
2Lok{) 



(C31) 



and the ratio of the angle deflection to the surface force 
is 



F, 



3C+ 
2Loko 



3L2S- 
Llko 



(C32) 



Using these, the effective spring constant, which is de- 
fined as the ratio of the vertical force to the vertical con- 
tact point deflection, is 



kcffifj.) = — 

Zc 

F, 



(C33) 



CoCi + SqL29l 



ko 



Co 



C4 



2Lo 



3C+ 
2Lo 



When one speaks of 'the' effective spring constant, one 
means out of contact, fccff = fccff(O), (which means 
C+ — Co, S- — So)- One can see that for the case 
L2 — Q; which means that Zc = CoCl, this is the same 
as Eq. (|C26p . These expression agree with those given in 
§S|B9]and[BT0l 

These three equations are sufficient to obtain every- 
thing that is needed from an atomic force measurement 
of surface forces that has been analysed to give the an- 
gular deflection. That is, they give Fz{9l), and hence 
Zc{Fz), and hence h = Zp — Zc{Fz) + zo- This assumes 
that the the angular deflection and the zero of separation 
have been obtained as described in JJB] It also assumes 
that any variation in distributed forces on the cantilever 
can be neglected, which means that these are apparent 
rather than real quantities. 



4. Drag Force on Cantilever 

In the atomic force microscope one measures a drag 
force at large separations 



-Fdrag = — 67r77Ldrag2^p- 



(C34) 
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This defines the drag parameter Ldrag- The drag length 
is an effective quantity that accounts for the cantilever 
shape, the fact that the drive velocity is not orthogonal 
to the cantilever axis, (if one wanted to, one could instead 

define i^drag = -BTrryLdragC'oip, with Zdrag = -'^drag/Co), 
and, most importantly, the fact that the distributed drag 
force is interpreted as if it were a point force on the con- 
tact point. 

The drag force per unit length may be taken to be 



/(ct) = -Cdrag[ipCo + Cicr)]. 



(C35) 



This simply says that the local drag force is proportional 
to the local velocity of the cantilever through the liquid. 
This drag force is orthogonal to the axis of the unde- 
flected cantilever, which was how the distributed force 
was defined. 

At large separations C.i'^) ~ and this force is uniform. 
In this case, Eq. (|C28p gives the deflection angle as 



2k'n ■ 



(C36) 



The end force and moment due to a point force, Fz = 
-Fdrag, and no friction, Fy = 0, is Fl = Co^drag and 
Ml = i2<5'o-Fdrag, wMch gives a deflection angle 



3Co 3L2S0 



2Loko L^ka 



F 



drag- 
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Equating these two expressions for the deflection angle 
gives 



Cdr 



2fco 
ZpC'o 



3Co , 3L2S0 



2Loko 



L^ko 



127r?7£drag 

Co 



3Co 
2Lo 
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3L2SQ 



^0 
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The leading term is identical to that used previously,^ 
(Previously the subordinate term due to tilt and torque 
were not accounted for.) 



5. Comparison of Theory and Measurement 

Next an algorithm will be given for computing the de- 
flection of the cantilever due to the drainage and van der 
Waals forces on the colloid probe and the drag force on 
the cantilever. This algorithm gives the actual vertical 
deflection Zc, the actual separation h, and the actual sur- 
face force Fz, as well as the actual angular deflection 6^. 
The atomic force microscope measures 9l, and assumes 
that this angular deflection is due solely to a point force 
Fz, and uses the point force equations to convert the 
angular deflection to an apparent vertical deflection, an 
apparent separation, and an apparent point force. This 
procedure would be accurate if the distributed drag force 
on the cantilever were constant, (because then the dis- 
tributed force contribution to the deflection is just a con- 
stant that is removed in the base line treatment of the 



data) . However to the extent that the drag force is vary- 
ing with cantilever deflection, the apparent quantities are 
not equal to the actual quantities. Therefore, in compar- 
ing theory with measurement one has to convert an exact 
angular deflection given by theory to an apparent verti- 
cal deflection, an apparent separation, and an apparent 
point force, using the same conversion factors that one 
uses to analyze the atomic force microscope raw data. 
(Out of interest, one can always compare the actual and 
the apparent calculated quantities to see how big an effect 
variable drag has.) In the text above, it is the apparent 
theoretical and experimental quantities that are plotted. 



6. Numerical Computation for Variable Drag 

The shape of the cantilever (^{a) under the influence 
of the point and distributed forces was given above as 
Eq. (jC23[) . The shape is a function Fz (via the fric- 
tion force Fy, if present, the end force F^, and the end 
moment Ml), which gives it a direct dependence on h 
and h, and a functional of C(c) (via f{cr))- This func- 
tion will be denoted S{a;Fz, [C]), or more simply S{a). 
(This is necessary because in the numerical solution of 
the equations of motion, one has to distinguish and to 
equate the shape that evolves from the previous shape, 
C(i-f At), and the shape that satisfles the elasticity equa- 
tion, S{a{t + IS.t))-) The deflection of the end for a given 
shape is = C(-^o) = S{Lq). Similarly the deflection 
angle of the end is 9l — diS'(cr)/dcr|^^^^. The separa- 
tion is h = Zp + Zc, and the rate of change of separation 
is h = Zp + Zc, where Zc is the function of Cl and 9l 
given above, Eq. (IC18I) . For convenience the constant of 
separation has been chosen as zero. 

The problem to be solved is: given the drive trajec- 
tory Zp{t), at each instant calculate the shape S{(T;t) 
and hence CL{t), dL{t), h{t) etc. The axial coordinate 
can be discretised, ai = zAo- and the shape can be writ- 
ten simply as a vector C,, with elements Q = C(o'i)- Also 
i e [0, L], (Jl = Lq and Cl = C,{vl) = C{La), consistent 
with the notation introduced above. 

Suppose the shape C, and its rate of change C, is known 
at time t. It is required to obtain these at time t' = i-l-Aj. 
One has three equations to solve: 



c = c 



A*C + :^A?C, 



C + AtC, 

^{F'zX')- 



(C39) 



Clearly then one has a single unknown vector, namely the 
acceleration C,, and one can insert the flrst two equations 
into the third, which can then be written as a function 
of the acceleration, C'(C) = S(C)- (Both sides of this 
depend upon the known ^ and C at the preceding time 
step.) Equivalently then one has to flnd the zero of the 
vector function 



G(0 = C'(C)-s(C). 



(C40) 
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In order to develop a stable iteration procedure to solve 
this, one converts it to a fixed point problem and uses 
essentially Newton's method. That is, define 

H(O^C-(^) 'g(C), (C41) 

where the Jacobean matrix has elements <^ dG /dC> ~ 

dGi/dC,j. Clearly the correct acceleration C,, which satis- 
fies G(C) = is a fixed point of this, H(^) = ^. Hence 
one can set up an iteration procedure for the acceleration, 

^("+1) =H(C^")). (C42) 

The procedure has been designed to be quadratically con- 
vergent by making the Jacobean of H vanish at the fixed 
point. To avoid matrix inversion, this is more conve- 
niently written 



'(n+l) 



g("). 



(C43) 



5C dC 

Care has to be taken with the matrix multiplications, 
since these should be weighted according to the trape- 
zoidal quadrature rule. Gaussean elimination can be used 
to obtain from this the next iterate for the acceleration. 
For the first guess, the converged acceleration from the 
preceding time step can be used, C^^Kt') = C{t). 
The elements of the Jacobean matrix are explicitly 



with f;, ^ F,{t + At) 
has 



-M^„A,^. (C44) 
A Kronecker-(5 appears here. One 



-A^Cdrag/S 



1 + 5jO + 5: 
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dF^ 2B dF^ 
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Co + p.So, and 



dF, 



where the effective friction constant was given above. Fi- 
nally, 
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The final term follows because the deflection angle of the 
end is 9l = [Cl — Cl-i]/^^. The force derivatives were 
given in 21 

Typically L = 100 points were used to discretise the 
axial integrals. This many points markedly slows the 
Gaussean elimination. Therefore an orthonomal basis for 
a sub-space was constructed, w„, with n = 1,2,..., with 
elements Vni & polynomial in yi of degree n. The polyno- 
mials are chosen orthogonal on the interval [0, Lq] with 
trapezoidal weighting. The orthonomal basis was con- 
structed numerically. Typically four vectors were used 
for the basis. The various matrices and vectors were pro- 
jected onto the sub-space and the iteration for the accel- 
eration carried out. Typically 5 iterations were used at 
each time step. Tests showed the results did not change 
upon increasing the number of iterations, basis vectors, 
or grid points. The procedure was found to be quite sta- 
ble. This includes in contact, the turn around point, and 
the jump from contact. 



7. Simplified Algorithm for Constant Drag 

The algorithm simplifies considerably if one can neglect 
the variation in cantilever drag with deflection. In this 
case one only needs the vertical deflection Zc, the velocity, 
ic, and the acceleration Zc of the contact point. With 
again t' = i 4- At, the three equations to be solved are 



1 



z'^ = Zc -I- Atic -f ^^t-^c, 

~ + AtZc, 

Zc = F^{h\h')/kc«. 



(C49) 



RecaU F^ = Fvdw + ^drain- Bccausc h' — Zp + z'^, the 
right hand side of the last equation is a function of Zc. 
Hence the correct acceleration is the zero of 



G(ze) = z^(zc) - F,{h', h')/k,s. (C50) 
This is converted into a fixed point problem by defining 



L2So-fiL2Co, (C47) H{2<: 



(C51) 



Ze + AtZe + O.SApe - F, jh' , h' ) / k^g 

A? dF,{h',h')A} dF,{h',h')^ 



2 hsdh' 2 k^gdh' 

From this, the iteration procedure for the acceleration is 

z("+i) = i7(z(")). (C52) 

This is quadratically convergent and was found to be 
quite stable, including in the contact region where the 
van der Waals force is steeply repulsive. 



